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In this paper, we consider quadratic Maxwell invariant as a correction to the Maxwell theory 
and study thermodynamic behavior of the black holes in Einstein and Gauss-Bonnet gravities. We 
consider cosmological constant as a thermodynamic pressure to extend phase space. Next, we 
obtain critical values in case of variation of nonlinearity and Gauss-Bonnet parameters. Although 
the general thermodynamical behavior of the black hole solutions is the same as usual Van der Waals 
system, we show that in special case of the nonlinear electromagnetic field, there will be a turning 
point for the phase diagrams and usual Van der Waals is not observed. This theory of nonlinear 
electromagnetic field provides two critical horizon radii. We show that this unusual behavior of 
phase diagrams is due to existence of second critical horizon radius. It will be pointed out that the 
power of the gravity and nonlinearity of the matter field modify the critical values. We generalize the 
study by considering the effects of dimensionality on critical values and make comparisons between 
our models with their special sub classes. In addition, we examine the possibility of the existence 
of the reentrant phase transitions through two different methods. 


I. INTRODUCTION 

It is well known that black hole behaves as a thermodynamic system and it can be interpreted with a physical 
temperature and an entropy [ij. Finding the connection between the laws of black hole mechanics and the laws of 
ordinary thermodynamics is one of the remarkable achievements of theoretical physics during the last forty years. 
Therefore, it is very natural to study various thermodynamic aspects of black holes, such as thermal stability, phase 
transition and black hole evaporation. Studying phase transition in black holes would be one of the fascinating topics 
in this regard because this phenomena plays an important role in order to explore thermodynamic properties of various 
systems near the critical points. The first attempt to investigate the phase transition of the black holes has been done 
by Hut and Davies Q. 

In general, there are four different approaches for studying the phase transition of the black holes, theoretically. 
First, the cosmological constant. A, is considered to be pressure of the system, in which its related conjugate quantity 
will be volume. By this consideration, the critical behavior can be studied through phase diagrams. Second, studying 
phase transition in black holes by using the Clausius-Clapeyron-Ehrenfest’s equations Q- Considering the analogy 
between the thermodynamic state variables and various black hole parameters (U O Q and P e-)- — $), puts us in a 
position to write down the Ehrenfest’s equations for the black holes Q and study their phase transition. Third, an 
alternative approach to investigate the phase transition was suggested by Ruppeiner in 1979 @ in which proposed 
a geometrical way to study thermodynamical phase transitions. Eourth, investigating the phase transitions of black 
holes through the canonical ensemble by calculating heat capacity 

In this paper, we are going to investigate the phase transition of black hole solutions in the asymptotically AdS 
spacetime by considering the first mentioned approach (the cosmological constant as a pressure of the system) in both 
Einstein and Gauss-Bonnet (GB) gravities. In context of AdS/CET correspondence, it was proposed that variation of 
A corresponds to variation of the number of the colors on boundary of the Yang-Mills theory with chemical potential 
interpretation [J |1|. There are two main reasons to investigate the asymptotically AdS black holes. Eirst reason 
is that the AdS/GET correspondence attracts attentions to the physics of asymptotically AdS black holes in recent 
years; the main focus is on understanding strongly coupled thermal field theories living on the AdS boundary. Even 
from a bulk perspective such black holes and their thermodynamics which exhibits various phase transitions, are quite 
interesting. Second, the behavior of the black hole phase transition in the asymptotically AdS spacetime is like the 
Van der Waals liquid/gas [I,®. and also, is different from those of black holes in the flat space EH- Studying 
thermodynamic behavior of black holes in an asymptotically AdS spacetime has been done first by Hawking and 
Page in 1983 [l^. After that, the critical behaviors of the black holes by including the cosmological constant as a 
thermodynamic pressure have been investigated in In this approach, the black hole mass M is considered as 
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the Enthal py of the system. Studying phase transition of black holes with Einstein gravity has been done in many 
literatures |l4j. In addition, the critical behavior of charged AdS-GB black holes has been investigated in [T5l - [T7l| . 

One the other hand, in electrodynamic point of view, the self-energy of a point-like charge has a divergency at the 
origin. In order to remove this singularity, Born and Infeld introduced an interesting kind of nonlinear electrodynamics 
(NED) in 1934 [l^. Coupling of NED with the gravity was first done by Hoffmann [I^. The effects of Born-Infeld 
NED coupled to the gravitational field have been studied in various contexts such as superconductors , wormholes 
[21I, [23 , static black holes and rotating black objects 24 1. In addition, another motivation for considering Born- 
Infeld NED comes from the fact that it naturally arises in the low-energy limit of the heterotic string theory [ 2 ^. 
Recently, two different Born-Infeld types of NED have been introduced by Soleng [2^ and Hendi 27 1. The Soleng 


Lagrangian has a logarithmic form and, like Born-Infeld theory, removes divergency of the electric field while the 
Lagrangian proposed by Hendi has an exponential form and does not cancel the divergency of the electric field but 
its singularity is much weaker than that in Maxwell theory. Investigation of black object solutions coupled to these 
two nonlinear fields has been done in [ 2 ^. The Lagrangian of mentioned Born-Infeld type nonlinear theories, for 
weak nonlinearity, can be written with the following form 


C{T) = - + /3jr2 + o (/ 32 ) , 


( 1 ) 


where /3 is proportional to the inverse value of square nonlinearity parameter in Born-Infeld-type theories, so it gets 
just positive values. In Eq. ([1]), is the Maxwell invariant, = d^Ay — d^A^j^ is the electromagnetic field 

tensor and A^ is the gauge potential. In addition, /3 denotes nonlinearity parameter which is small. For P —> 0, C{F) 
reduces to the standard Maxwell Lagrangian, C-Maxweii {^) = as it should be. In this paper, we take into account 
the Eq. Ill) as a NED source coupled to the Einstein and GB gravities and investigate the effects of nonlinearity on 
the properties of the phase transition. 

It is worthwhile to mention the motivations for considering the NED Lagrangian and specially Lagrangian CD- 
Nonlinear field theories are of interest to different branches of mathematical physics because most physical systems 
are inherently nonlinear in the nature. The main reason to consider NED comes from the fact that these theories 
are considerably richer than the Maxwell field and in special case they reduce to the linear Maxwell theory. Various 
limitations of the Maxwell theory, such as description of the self-interaction of virtual electron-positron pairs and 
the radiation propagation inside specific materials [s^, motivate one to consider NED [ 2 ^. Besides, NED improves 
the basic concept of gravitational redshift and its dependency of any background magnetic field as compared to the 
well-established method introduced by standard general relativity. In addition, it was recently shown that NED 
objects can remove both of the big bang and black hole singularities [mI - [^ . Moreover, from astrophysical point of 
view, one finds that the effects of NED become indeed quite important in superstrongly magnetized compact objects, 
such as pulsars and particular neutron stars (also the so-called magnetars and strange quark magnetars) [s^. Also, 
since the gravitational redshift of magnetized compact objects is connected to the mass-radius relation of the objects, 
it is important to note that NED affects the mass-radius relation of the objects. It is worthwhile to mention that one 
can find regular black hole solutions of the Einstein field equations coupled to a suitable NED . In addition, 

an interesting property which is common to all the NED models is that these models satisfy the zeroth and first 
laws of black hole mechanics. The appropriate world-volume dynamics on a curved D3-brane may provide a plausible 
frame-work at Planck scale by incorporating the Einstein-NED. At this point, elimination of strong intrinsic curvature 
in the regime by the strong nonlinearity in the 17(1) gauge theory is remarkable [^. From the point of view of 
AdS/GFT correspondence in hydrodynamic models, it has been shown that, unlike gravitational correction, higher- 
derivative terms for abelian fields in the form of NED do not affect the ratio of shear viscosity to entropy density 
[ 3 ^ . Motivated by the recent results mentioned above and the fact that we accepted NED as a generalization of the 
Maxwell theory, it is natural to apply NED theories for charged objects such as black holes. 

Now, we focus on motivations of considering the nonlinear term of the electromagnetic field perturbatively. Although 
various theories of NED have been created with different primitive motivations, only for the weak nonlinearity (Eq. 
CD), they contain physical and experimental importances. As we know, using the Maxwell theory in various branches 
leads to near accurate or acceptable consequences. So, in transition from the Maxwell theory to NED, the logical 
decision is to consider the effects of weak nonlinearity variations, not strong ones. This means that, one can expect 
to obtain precise physical results with experimental agreements, provided one regards the nonlinearity as a correction 
to the Maxwell field. On the other hand, several reasonable papers have been published by considering Eq. CD as 
an effective Lagrangian of electrodynami cs 1^ . I37l - l43 |. Heisenberg and Euler have shown that quantum corrections 
lead to nonlinear properties of vacuum [29j. Also, it was proved that in the low energy limit of heterotic string 
theory, a quartic correction of the Maxwell field strength tensor appears [s^- So it is natural to consider Eq. CD 
as an effective and suitable Lagrangian of electrodynamics instead of the Maxwell one. Investigating the effects of 
nonlinearity parameter of Eq. CD coupled to the Einstein, GB and third order Lovelock gravities have been done 
in [H, OSlI and [i^, respectively. In this paper, we are dealing with stringy corrected electrodynamics, and 
therefore, obtained results are applicable in context of string theory. In order to separate the valid domains of the 
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stringy corrected and classical Born-Infeld theories, we will present a limiting point (vertical line in the fignres). 

Motivations for considering GB gravity can be found in literature. For example, we refer the reader to the interesting 
nontrivial causal structure of GB gravity with superluminal graviton modes |43l |. Finally, we should note that although 
most of thermodynamic works in black hole physics related to asymptotically AdS solutions, there were some attempts 
to investigate thermodynamical behavior of dS black holes [43 |. 

The outline of our paper is as follows. Section II is devoted to introduction to Einstein and GB black hole 
solutions and their conserved quantities. Next, we extend the phase space by considering cosmological constant as 
thermodynamic pressure and calculate critical values and then we plot diagrams for different cases. We give a detailed 
discussion regarding diagrams, their physical interpretations, and the effects of both nonlinear electromagnetic and 
gravitational parameters. We finish our paper with some closing remarks. 


II. FIELD EQUATIONS AND CONSERVED QUANTITIES 


In order to study phase transition of Einstein-Gauss-Bonnet gravity in presence of a generalized nonlinear electro¬ 
magnetic field, one can employ the following Lagrangian 

•^tot ~ ^EN 2K + aCGB+C{T), (2) 

where the Lagrangian of Einstein gravity is the Ricci scalar, Cen = and A is the negative cosmological constant. 
In third term of Eq. m, a is the GB coefficient with dimension (Length)^ and Cgb is the Lagrangian of GB gravity 
with following form 

Cgb = RabcdR'^^^'' - 4.RabR‘^^ + . (3) 

Using variational method, we obtain the following field equations 

Gf, + kgab + = \gabC{T) - (4) 


da {V^CeF‘^^) = 0, (5) 

where Gfj, is the Einstein tensor, = 2 {RacdeRt‘^''- 2RacbdR‘'‘^ - 2RacRb + R-Rab) - ll2CGBgab and Ljr = 

dC{T)ldT. 

Now, we are interested in studying topological black holes and their phase diagrams, therefore, we employ the 
following static metric 

dr"^ 

ds^ = -f{r)df + + r^dnl_j^, (6) 


in which 


dni 


n-1i-1 

^^1 + S n k = 1 

i=2 j=l 

n —1 2—1 

< dOf + sinh^ Oidd^ + sinh^ X) 0 sin^ djddf k = —1 

i=3 j=2 

n—1 

dQf fc = 0 


(7) 


with volume Un-i- 

We use Eq. m and mentioned metric to obtain radial electromagnetic field tensor as 


Ftr — 


-.n— 1 


7^3n—3 


o{p^). 


( 8 ) 


In order to find Einstein solutions one can use two methods: one, by putting a = 0 and using mentioned field 
equations. Second approach is obtaining GB metric function through use of field equations and series expanding it 
for small values of GB parameter. In order to give more specific details, we use the second approach. Therefore, in 
case of GB gravity, one can obtain metric function in form of [4lj 


2 

/(r) = k + ^{l - V'4'(r)) 


(9) 








with 
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'l'(r) = 1 + 


8 a' 


n{n — 1 ) 


A- 


i(n — l)r 

2rn 


nq 


2nq'^f] 


(n — 2 )r^ 


j,4n-4('3^ _ 4^ 


+ O (/32) 


( 10 ) 


where m is an integration constant that is related to mass and a' = {n — 2)(n — 3)a. It is evident that for case of 
small values of nonlinearity, the metric function will lead to the GB-Maxwell gravity. As for Einstein gravity, series 
expanding of GB metric function for small values of a' will lead to 


/ {r) = Jen - 




(n — 1) (3n — 4) 


-6 


/3+%^a' + 0(a'/3,a'^/32), 


( 11 ) 


where the metric function of Einstein-Maxwell gravity is 

fEN = k- 


m 


2q^ 


.(n- 1) 


(n — 1 ) (n — 2 ) r 


2n-4 • 


( 12 ) 


Next step is devoted to calculating conserved quantities. In general, for both Einstein and GB gravities one can 
find total mass of black hole in form of 1 ^, 


M = 


Wn-I (n - 1 ) m 
1071 


(13) 


It is notable that, although the form of total mass in GB and Einstein gravities seems to be the same, its value is 
different for Einstein and GB branches (fc ^ 0). In other words, the geometrical mass for GB gravity will be 


rriGB = k [ka + r\) r" - 


2r"A 
n (n — 1 ) 


2q^ 




{n — 1 ) (n — 2 ) r_|_ 


n-2 


(n - 1) (3n - 4) 


-+0(/3") 


(14) 


where r+ satisfies /(r = r+) = 0, and Eq. (TT^ reduces to geometrical mass for Einstein gravity in case of a' = 0. 

Previously, it was seen that obtained metric functions are representing black holes with essential singularity located 
at r = 0. Geometrical properties of the solutions were investigated in [sl, |4l| and it was shown that these solutions 
can be interpreted as asymptotically AdS black holes. Therefore, by using the definition of surface gravity and its 
relation with Hawking temperature we find temperature of these two black holes as [4T| 


T = 


kin - l)(n - 2)rX--^ (l + - 2 rt"-^A - 2 r^"- 2 g 2 + 4 ^ 4 ^ 

47r(n-l)rt"-^(l + ^) 


O 


(15) 


where in order to find temperature of Einstein gravity, it is sufficient to set a' — Q 

Due to the fact that solutions are asymptotically AdS, in order to find entropy of these two gravities, one can use 
Gibbs-Duhem relation. Therefore, we obtain the following relation [4l| 


5 = 


Vn-l 

4 



(n-3)4 ) 


jn — 1 


(16) 


where in order to find entropy related to Einstein gravity, one should set GB parameter to zero [ 3 ^ . 

Generally, as one can see, the topological structure of spacetime modifies the amount of contribution of GB pa¬ 
rameter, for the GB gravity. Although the presence of GB gravity in flat case is evident in metric function, regarding 
to conserved and thermodynamic quantities no contribution of GB gravity was seen and obtained values are same 
as Einstein gravity. Considering this fact and equation of Gibbs free energy, the thermodynamic behavior and phase 
diagrams of GB and Einstein gravities in case of flat horizon are same. In addition, it was shown that [H, |4l| obtained 
conserved and thermodynamic quantities satisfy the first law of thermodynamics with the following form 


dM = TdS -b A>dQ, 


(17) 


where Q = ^ and 4> = J Ftrdr\r=r+- 
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III. EXTENDED PHASE SPACE AND PHASE DIAGRAMS 


In order to investigate the phase structure of the solutions, we employ the approach in which the cosmological 
constant is a thermodynamic variable corresponding to thermodynamical pressure with the following relation 


P = - 


Stt 


(18) 


This consideration could be justified due to the fact that in quantum context, fundamental fixed parameters could 
vary and they are not fixed. In the absence of cosmological constant, a sourceless solution of the Einstein’s equation 
is Minkowski spacetime where the isometry transformations are governed by Poincare group. In the presence of 
cosmological constant, Minkowski is no longer a valid solution and it is replaced by the (anti-)de Sitter spacetime with 
(anti-)de Sitter group description. Considering that we are now employing the (anti-)de Sitter group for describing 
kinematics, the ordinary notions of energy and momentum, as well as the relationship between them and the causal 
structure of spacetime will be modified [^. Therefore, we expect to see the effects of this distortion on spacetime of 
black holes and its corresponding thermodynamical values which is evident from calculated thermodynamical values. 
As one can see the conjugating thermodynamic variable to this assumption (cosmological constant as pressure) will 
be volume where in literature the derived volume for different types of black holes are same as that for the topology 
of the spacetime [3, [lH- order to calculate the volume of these thermodynamical systems, we use the following 
relation 


V = 



(19) 


In addition, it was proven that the Smarr formula should be extended to Lovelock gravity as well as nonlinear 
theories of electrodynamics [3, [3) Sll ■ Scaling argument was used to derive an extension of the first law and its 
related modified Smarr relation that includes variations in the cosmological constant. Lovelock coefficient, and also 
the nonlinearity parameter [3)[3Ell- In our case, perturbative Lovelock gravity with NED, M can be a function of 
entropy, pressure, charge. Lovelock parameter, nonlinearity coupling coefficient. Regarding the previous section, we 
find that those thermodynamic quantities satisfy the following differential form. 


dM = TdS + ^dQ + VdP + Ada + Bdfd, 


( 20 ) 


where 


UM\ 


S,Q,P,/3 


{n — l)fc^r(|_ 
IGtt 


n—4 


B = 


J S,Q,p,a 4(3n - 4)7rr: 


,3n-4 ■ 


Moreover, scaling argument helps us to obtain the generalized Smarr relation for our black hole solutions in the 
extended phase space 


77 — 1 

M = - -TS + 4>dQ 

n — 2 



2(n-l) 

(n-2)(n-3)'^“ + 



( 21 ) 


We should note that the obtained relations is valid for perturbative Lovelock gravity with NED. Hereafter, we treat 
the cosmological constant as a thermodynamical variable, while Gauss-Bonnet coefficient and nonlinearity parameter 
of NED as two constants. With doing so the total finite mass of the black hole will play the role of Enthalpy and the 
corresponding Gibbs free energy will be in form of 


G = H-TS = M-TS. 


( 22 ) 


The obtained volume for our considered cases is 


V = 


UJn-ir+’" 

n 


(23) 


which is consistent with topological structure of spherical symmetric spacetime. This result is consistent with what 
was derived previously [14l . Il5l | and shows the fact that although considering GB gravity modifies the metric function 
and some conserved quantities of the black hole, it does not change the volume of the black hole. In other words, the 
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volume of the black hole is solely dependant on the cosmological constant (pressure). Due to relation between volume 
and radius of the black hole, we use horizon radius (specific volume) in order to investigate the critical behavior of 
these systems 

Next step will be calculating critical values. In order to do so, we use the method in which critical values are 
obtained through the use of P — r-|_ diagrams. Since the critical point is an inflection point on the critical isotherm 
P — r+ diagram, we use the following relations to obtain the proper equations for critical quantities 


dP 

dr+ 




(24) 


It will be constructive to give a short description regarding different phase diagrams and the information they 
contain before presenting tables and phase diagrams. G — T diagrams are representing energy level of different states 
that phase transition takes place between them. The characteristic swallow tail that is seen in these diagrams shows 
the process that we know as phase transition. It also gives interesting information regarding temperature of critical 
points. For T —r_|_ diagrams, it contains information regarding critical temperature and horizon radius in which phase 
transition takes place. Also, it gives some insight about single state regions which in our case is small/large black 
holes. It also helps us to understand the effects of different parameters on critical temperature and horizon radius, 
and whether by changing value of a parameter, system needs more or less energy in order to have phase transition. 
If one is interested in studying conductor/superconductor transition that these nonlinear electromagnetic fields are 
representing, studying these diagrams will give more information regarding to conductivity and superconductivity 
regions. Studying P — diagrams gives us information regarding the behavior of pressure as a function of horizon 
radius, and critical pressure and horizon radius of phase transition. Finally, using the fact that the free energy, 
temperature, and the pressure of the system are constant during the phase transition, one can plot the coexistence 
curve of two phases. One of the reasons for studying these diagrams is the similarity between phase structure of 
black holes and the Van der Waals thermodynamical systems. Here, we have used the geometric units and investigate 
thermodynamic behavior and critical point of the solutions, qualitatively. 

Using Eq. one can find Tc in one of the equations and replace it in other equation which leads to the following 
relations for calculating critical horizon radius 

{ fc (n — 2) + 8 (4n — 5) — 2 (2n — 3) = 0, Einstein 

(12n — 48) a'^k^ — 12a'k‘^r^ + {4 (4n — 5) — (2n — 3) 2q^ ■ (25) 

GB 

+ [{4 (4n — 7) — (2n — 5) 129 ^ 0 ;' + (n — 2) r^] k = 0, 

As one can see, due to complexity of the obtained relation for critical horizon radius, it is not possible to find critical 
horizon radius analytically. Therefore, we employ the numerical method in order to calculate critical quantities and 
study the effects of variation of parameters in case of spherical horizon {k — 1). The numerical calculations show 
that for cases of flat {k = 0) and hyperbolic (fc = —I) horizons, the critical pressure and critical temperature for both 
Einstein and GB gravities are negative, so like GB-Maxwell black holes [l^ , the phase transition does not take place. 
Here, we present various tables in order to study the effects of different parameters on critical values. Next, by using 
the information of these tables, we plot P — r+, T — and G — T diagrams for Einstein and GB gravities in the 
presence of nonlinear corrected Maxwell field (Figs. [T]~ ITUl) . It is notable that following results for critical pressure 
and temperature are obtained by using larger critical horizon radius. In order for higher orders of corrections to be 
small enough and do not acquire values higher than Maxwell term, we have plotted a vertical line which represents 
the limit for different cases. 

It is notable to mention that in order to have a well-defined vacuum solution with m = q = 0, the pressure P has 
to satisfy the following constraint 


^ 647ra'P 

n[n — 1) 

which puts a large bound for the pressure as maximal pressure 

n(n — 1) 


P < P — 

J max — 


647rQ!' 


(26) 


(27) 


It means that only for sufficiently small pressures, the solution Eq. m possesses an asymptotic AdS region. 
However, in this paper, we choose suitable parameters for pressure to be smaller than maximal pressure everywhere. 
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FIG. 1: Maxwell solutions for Einstein gravity: P — r+ for T — Tc (Left), T — r+ for P = Pc (Middle) and G — T for 
P = 0.5Pc (Right) diagrams for fc = l, n = 3 (continuous line), n = 4 (dotted line) and n = 5 (dashed line). 
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FIG. 2: Maxwell solutions for GB gravity: P — r+ for T = Tc (Left), T — r+ for P = Pc (Middle) and G — T for P = O.SPc 
(Right) diagrams for k = 1, q — 1, a' = 10“^, n = 4 (continuous line), n = 5 (dotted line) and n = 6 (dashed line). 




FIG. 3: Maxwell solutions for Einstein and GB gravities: P — r+ for T = Tc (Left), T — r+ for P = Pc (Middle) and 
G — T for P = O.SPc (Right) diagrams for k = 1, n = 4, q — 1, a' = 0.1 for Einstein (continuous line) and GB (dashed line) 
gravities. 















FIG. 4: Einstein solutions: P — r+ (Left), T — r+ (Middle) and G — T (Right) diagrams for k — 1, n — 4, q = l and 
/3 = 0.045. 

P — r+ diagram: T = O.STc (continuous line), T = Tc (dotted line) and T = 1.2Tc (dashed line). 

T — r+ diagram: P = 0.5Pc (continuous line), P = Pc (dotted line) and P = 1.5Pc (dashed line). 

G — T diagram: P — 0.5Pc (continuous line), P = Pc (dotted line) and P = 1.5Pc (dashed line). 




T 


FIG. 5: Einstein solutions: P — r+ (Left), T — r+ (Middle) and G — T (Right) diagrams for fc = l, n = 4, and j3 = 0.07. 

P — r+ diagram: T = O.STc (continuous line), T = Tc (dotted line) and T = 1.2Tc (dashed line). 

T — r+ diagram: P — 0.5Pc (continuous line), P = Pc (dotted line) and P = 1.5Pc (dashed line). 

G — T diagram: P = 0.5Pc (continuous line), P — Pc (dotted line) and P = 1.5Pc (dashed line). 


f3 G Tc Pc ^ 


f3 rc Tc Pc ^ 

0 2.44948 0.04330 0.00331 0.18750 


0 1.49534 0.17029 0.03558 0.31250 

0 3120Q 

0.00100 0.04331 0.00332 0.18748 

2.44917 


0.00100 0.17030 0.03559 0.31247 

1.49505 

0.05000 0.04343 0.00333 0.18569 

2.43312 


0.05000 0.17110 0.03596 0.31093 

1.47940 

1 03212 

0.10000 0.04356 0.00336 0.18635 

2.41568 


1 01944 

0.10000 0.17201 0.03639 0.30890 

1.46018 

1 1 

0.15000 0.04370 0.00338 0.18569 

2.39698 


1 11349 

0.15000 0.17310 0.03690 0.30603 

1.43543 


Table 1 (left): Einstein gravity for q = I and n = 3. Table 2 (right): Einstein gravity for g = I and n = 4. 
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FIG. 6: Einstein solutions: P — r+ for T — Tc (Left), T — r+ for P — Pc (Middle) and G — T for P = 0.5Pc (Right) diagrams 
for fc = 1, n = 4, g = 1, /3 = 0.045 (continuous line) and /3 = 0.07 (dashed line). 





FIG. 7: GB solutions: P — r+ for T — Tc (Left), T — r+ for P = Pc (Middle) and G — T for P = 0.5Pc (Right) diagrams for 
k = 1, n = 4, q — 1, a' — 10“^, and /3 = 0 (continuous line), p = 0.03 (dotted line) and P = 0.06 (dashed line). 


/3 re Tc Pc ^ 


/3 re Tc Pc ^ 

0 1.29271 0.31658 0.10714 0.43750 


0 1.49550 0.17024 0.03557 0.31248 

0 55202 

0.00100 0.31661 0.10716 0.43746 

1.29247 


0.01000 0.17040 0.03564 0.31219 

1.49252 

0 91471 

0.05000 0.31804 0.10820 0.43532 

1.27952 


0.03000 0.17071 0.03579 0.31158 

1.48626 

1 01396 

0.10000 0.31978 0.10949 0.43221 

1.26233 


0.05000 0.17105 0.03594 0.31091 

1.47957 

1 09170 

0.15000 0.32202 0.11117 0.42668 

1.23597 


0.07000 0.17140 0.03610 0.31017 

1.47237 


Table 3 (left): Einstein gravity for g = 1 and n = 5. Table 4 (right): GB gravity for g = 1, a = 10 and n = 4. 
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FIG. 8: GB solutions: P — r+ (Left), T — r+ (Middle) and G — T (Right) diagrams for fc = l, n = 4, (/=! and j3 = 0.07. 
P — r+ diagram: T — Tc and (x = 0.05 (continuous line), a' = 0.1 (dotted line) and (x = 0.5 (dashed line). 

T — r+ diagram: P = Pc and a' = 0.05 (continuous line), a' = 0.1 (dotted line) and a' = 0.5 (dashed line). 

G — T diagram: P = 0.5Pc and a' — 0.05 (continuous line), a' = 0.1 (dotted line) and a' = 0.5 (dashed line). 





FIG. 9: Einstein solutions: P — r+ for T — Tc (Left), T — r+ for P = Pc (Middle) and G — T for P = 0.5Pc (Right) diagrams 
for fc = 1, g = 1, /3 = 0.001, n = 7 (continuous line), n = 8 (dotted line) and n = 9 (dashed line). 


a 

Tc 

7^ P Pc'^c 

c ^ c 

0.01000 

0.94914 

1.48981 

0.16688 0.03458 0.30873 

0.05000 

0.94608 

1.55579 

0.15157 0.02955 0.30333 

0.10000 

0.94374 

1.63116 

0.13707 0.02499 0.29747 

0.50000 

0.93836 

2.12855 

0.08560 0.01089 0.27103 

1.00000 

0.89005 

2.65829 

0.06346 0.00619 0.25942 


/3 

Tc 

rp p PcVc 

c ^ c Y 

0 

1.49550 

0.17024 0.03557 0.31248 

0.01000 

0.67158 

1.49252 

0.17040 0.03564 0.31219 

0.03000 

0.81326 

1.48626 

0.17071 0.03579 0.31158 

0.05000 

0.89202 

1.47957 

0.17105 0.03594 0.31091 

0.07000 

0.95018 

1.47237 

0.17140 0.03610 0.31017 


Table 5 (left): GB gravity for g = 1, /3 = 0.07 and n = 4. Table 6 (right): GB gravity for g = 1, a = 10 and n = 5 
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FIG. 10: GB solutions: P — r+ for T = Tc (Left), T — r+ for P = Pc (Middle) and G — T for P = 0.5Pc (Right) diagrams 
for fc = 1, g = 1, /? = 0.06, a' = 0.5, n = 4 (continuous line), n = 5 (dotted line) and n = 6 (dashed line). 



FIG. 11: Goexistence line for Einstein (Left) and GB (Right) gravities: P — T for A: = 1, n = 4 and g = 1. 
Left diagram: /3 = 0.045 (continuous line) and /3 = 0.07 (dotted line). 

Right diagram: a' = 10““* and /3 = 0 (continuous line) and /3 = 0.06 (dotted line). 


a 

Tc 

rp p Pc^c 

^ C ^ C T 

J- r. 

0.01000 

0.95884 

1.28530 

0.30889 0.10367 0.43139 

0.05000 

0.95598 

1.32895 

0.27681 0.08766 0.42089 

0.10000 

0.95414 

1.37662 

0.24745 0.07358 0.40937 

0.50000 

0.95469 

1.65361 

0.15000 0.03196 0.35241 

1.00000 

0.96227 

1.91323 

0.11049 0.01837 0.31815 


Table (7): GB gravity for g = 1, /3 = 0.07 and n = 5. 

It is a well-known fact that inclusion of higher order polynomial terms (e.g. rlf., d G N) in the Van der Waals 
equation of state changes the location of critical point in thermodynamic space {P,V,T), but does not change the 
universality class and therefore leads to the same exponents as the Van der Waals fluid (see [H, [l^ for more details). 
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FIG. 12: Coexistence line for Einstein (Left) and GB (Right) gravities: P — T for fc = 1, n = 4 and q = 1. 
Left diagram: /3 = 0.001 and n = 7 (continuous line), n = 8 (dotted line) and n = 9 (dashed line). 

Right diagram: (3 = 0.06, a' = 0.5 and n = 4 (continuous line), n = 5 (dotted line) and n = 6 (dashed line). 



FIG. 13: Coexistence line for GB gravity: P — T for A: = 1, n = 4, g = 1, /3 = 0.07 and a' = 0.05 (continuous line), a' = 0.1 
(dotted line) and a' = 0.5 (dashed line). 


This is different from the case of pure Lovelock gravity solutions [1^, which the equation of state is no longer a 
polynomial form. 


IV. DISCUSSION ON THE RESULTS OF DIAGRAMS 

In order to study the behavior of phase transition for these black holes in more details, we have plotted P — r_|_, 
T — r+ and G — T diagrams. For having better insight regarding the effects of correction on critical behavior of the 
system, we have also plotted some diagrams for the case of /3 = 0 which is the Maxwell theory. 

like Van der Waals system, the usual characteristic swallow tail is seen in black holes only in the presence of linear 
Maxwell field (right panels of Figs. [T]and[2]). In other words, considering Maxwell electromagnetic field leads to a 
Van der Waals like behavior and usual phase transition. The existence and usual behavior of phase transition are 
also evident from studying P — r+ and T — r+ diagrams. In case of the absence of nonlinearity parameter, for P — r+ 
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(left panels of Figs. [T]and[2]), pressure is a decreasing function of horizon radius but for a range of horizon radius, it 
is an increasing function of it which is not a physical behavior. Before and after this region there are two values of 
horizon with same pressure. The phase transition takes place between these two points which in case of black holes, 
it is small/large black hole phase transition. It is crucial to mention that, this region is only seen for case of T < Tc 
whereas for T > Tc pressure is only a decreasing function of horizon radius. On the other hand, for /3 = 0 in case 
of T — r+ diagrams (middle panels of Figs. [T]and[2]), if P = Pc the temperature is an increasing function of horizon 
radius and for a region of horizon radius, temperature is fixed. This place is the region where phase transition takes 
place and known as subcritical isobar. 

As for the effects of dimensions on the critical behavior of the system in the presence of linear Maxwell field Figs. 

□ and [5] are plotted. As one can see, the swallow tail (right panels of Figs. [I]and[2]), critical pressure (left panels of 
Figs. [Hand[5]) and temperature (middle panels of Figs. |T]and|5]) are increasing functions of dimension while critical 
horizon radius and subcritical isobars are decreasing functions of it. 

In addition, for the case of absence of nonlinearity parameter, we compare critical behavior of these two gravities 
with each other (Fig. |3|). As one can see considering GB gravity leads to increasing the size of swallow tail (right 
panel of Fig. El). The place of swallow tail is also shifted to lower values of temperature. On the other hand, pressure 
(left panel of Fig. [3]) and temperature (middle panel of Fig. [3]) of the critical point are greater in Einstein gravity 
whereas the length of subcritical isobar and the critical horizon radius are greater in GB gravity. Also, the needed 
energy for the phase transition in Einstein gravity is more than GB case. 

Next, we are considering /3 ^ 0 and plot one set of graphs for Einstein gravity for variation of nonlinearity parameter 
(Figs. 131131) and Figs. [7] and [5] for variation of nonlinearity and GB parameters for GB gravity. Also, we plot Figs. 
M and uni to investigate the effects of dimensions on the solutions of Einstein and GB gravities, respectively. It is 
notable that in these figures physical solutions exist only after the vertical line. This means that before these vertical 
lines the second term in Eq. ([5]) is not small enough (with respect to Maxwell term). The vertical lines in P — r+ 
(T —r+) and G — T diagrams are interpreted as the minimum value of the authorized horizon radius and temperature, 
respectively. In comparing figures there are more than one vertical line. Because of overlapping of these vertical lines 
with the vertical line of G — T diagrams, these lines are not presented. 

In case of variation of nonlinearity parameter same behavior is observed for Einstein and GB gravities. For this case 
the following results are obtained. Interestingly, contrary to Maxwell theory, in this case (/3 ^0), the Van der Waals 
like behavior is not preserved. The plotted graphs for Gibbs free energy versus temperature show the existence of a 
phase transition, a turning point. In other words, the characteristic swallow tail of phase transition in this nonlinear 
theory is modified and its shape is different from the usual thermodynamical systems. For small values of nonlinearity 
parameter (right panel of Fig. 11, the usual swallow tail is observed with a turning point and a minimum temperature 
(vertical line) which are located before swallow tail. It is evident that the minimum energy and the distance between 
turning point and swallow tail, are decreasing functions of /3 whereas the minimum temperature is an increasing 
function(see Figs. ElandlT]). It is worthwhile to mention that minimum temperature and critical temperature are two 
different quantities. 

As for the P — r_|_ , it is evident that the related graphs are modified like G — T diagrams. First, pressure is an 
increasing function of horizon radius, then after a turning point, it changes into being a decreasing function of r_|_. 
In this case, a part of the graphs shows the usual behavior of phase transition whereas there is another part which 
is irregular (left panels of Figs. |4]and|5l). In case of T = Tc, one can find two horizon radii for critical pressure. In 
this case the critical horizon radius is a decreasing function of nonlinearity parameter whereas the turning point and 
critical pressure are increasing functions of it (Tables 1 — 4 and 6). 

Moreover, in studying T — r+, same abnormal behavior is observed (middle panels of Figs. [4] and [5]). Usually, we 
are expecting temperature to be a decreasing function of horizon radius except in place of phase transition in which 
temperature is fixed and horizon radius increases. This region is know as subcritical isobars. But in this case, first 
temperature is a decreasing function of then it becomes an increasing function of it. As one can see in case of 
P = Pc the subcritical isobar is observed but another value of horizon radius exists which has the same temperature 
as subcritical isobar. This may show that in this place phase transition takes place. The critical temperature is an 
increasing function of nonlinearity parameter. 

Comparing the variational effects of GB and nonlinearity parameters, we find that they have opposite effects (Figs. 

□ andlS]). Thus, we leave out discussions of GB parameter for reasons of economy. 

It is worthwhile to mention a few characteristic behavior of graphs. As one can see, in case of small values of 
nonlinearity parameter, the distance between the critical point and turning point is large (right panels of Figs. Ejand 
0. Similarly, in case of P — r+ (T — r_|_) the region of r_|_ in which pressure being increasing (decreasing) function of 
r_|_ is large too. As nonlinearity parameter increases, the distance between these points decreases in G — T diagrams, 
and interestingly, in case of P —(T—r_|_) the region of r+, distance between the critical horizon radius, and turning 
point decrease (Figs. |3|and|71). In the end, we mention that PcTc/Tc is a decreasing function of nonlinearity and GB 
parameters. 
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FIG. 14: P versus r+ diagrams for q = 1 and fc = 1. 

left panel: a = 0 and /? = 0 (bold continues line), P = 0.03558 (continues line), /3 = 0.05 (bold dotted line), P = 0.03596 
(dotted line), /3 = 0.15 (bold dashed line), P = 0.03690 (dashed line). 

middle panel: a. = 10”"^ and /3 = 0 (bold continues line), P = 0.03557 (continues line), = 0.03 (bold dotted line), P = 0.03579 
(dotted line), /3 = 0.07 (bold dashed line), P = 0.03610 (dashed line). 

right panel: j3 = 0.07 and a = 0.1 (bold continues line), P = 0.02499 (continues line), a = 0.5 (bold dotted line), P = 0.01089 
(dotted line), a — 1 (bold dashed line), P = 0.00619 (dashed line). 


Next, we have considered the effects of dimensions on the critical behavior. As one can see, the temperature of 
swallow tail formation, energy gap between two states and minimum of temperature (turning point), critical horizon 
(middle panels of Figs. and fTOl). pressure (left panels of Figs. 1^ and ITOt and Pcfc/Tc are increasing functions of 
dimensions. An abnormal behavior for 4-dimension is observed in G — T diagrams for Einstein gravity. This behavior 
is due to power of in the last term of Gibbs free energy. The same abnormality could be obtained for the case of 
GB gravity in 5-dimension which is due to structure of Gibbs free energy, temperature and pressure. In other words, 
there are terms in these equations that vanish in case of n = 4. 

Finally, we have plotted the coexistence line in which along this curve, small and large black holes have alike 
temperature and pressure (Figs. [HI [121 and [13]). Gritical points are located at the end of the coexistence line where 
above these points the phase transition does not occur. In addition. Fig. [TT] indicates that nonlinearity parameter 
does not significantly affect the coexistence line. For both Einstein and GB gravities as dimension increases the critical 
pressure and critical temperature increase too (Fig. 1121) . In case of GB parameter, the critical pressure and critical 
temperature are decreasing functions of GB parameter (Fig. IT^and table 5). 


V. PHASE TRANSITION POINTS THROUGH HEAT CAPACITY 

Since we have observed an abnormal behavior in plotted phase diagrams, it will be worthwhile to test the existence 
of reentrant phase transitions through another method. In Ref. a new method for studying critical behaviors was 
introduced. The method is based on obtaining a relation for thermodynamical pressure by using the denominator of the 
heat capacity. In other words, by replacing the cosmological constant with its corresponding thermodynamical pressure 
and solving the denominator of the heat capacity with respect to this pressure a relation is obtained. This relation is 
different from the pressure which is obtainable from Eq. m- The maximums of this relation are representing places 
in which phase transitions take place. Therefore, the maximums in plotting a P — r+ diagram for this relation will 
give us critical pressure and horizon radius. This method was employed to obtain critical pressure in several papers 
which has proven to be an successful one [1^, . 

The heat capacity is obtained by 



-1 

Q 


(28) 


Using Eqs. (USD and ([ISD and replacing cosmological constant with its corresponding relation with pressure and 
solving its denominator with respect to pressure will lead to a relation for pressure which for economical reasons we 
will not bring it. We will present the results of this relation by considering mentioned values in different tables in 
following diagrams (Fig. [HI) . 
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First of all, in absence of the nonlinearity, the plotted diagrams shows existence of only one maximum. The pressure 
and horizon radius of this maximum is exactly located where phase transition takes place. Interestingly, by adding 
the nonlinearity to the system, two extrema are observed: a minimum and a maximum. The place of this maximum 
is exactly where phase transition occurs for mentioned vales for different parameters. But as one can see, no other 
maximum exists which indicates that, at least no other second order phase transition is observed. In other words, 
although for specific critical pressure in these diagrams two critical horizons are observed, the type of these critical 
points are not the same. Therefore, one can conclude that no reentrant phase transitions happens which is consistent 
with what was observed in coexistence diagrams. 


VI. CONCLUSIONS 

In this paper we have considered a quadratic Maxwell invariant as a correction term to the Maxwell Lagrangian. 
We studied the thermodynamic behavior of these solutions in Einstein and GB gravities. We considered cosmological 
constant as thermodynamic pressure and related conjugated quantity as volume of the black hole. By doing so, the 
interpretation of the mass as internal energy was changed into Enthalpy of the system. Therefore, not only the 
interpretation of the mass of the black hole was changed, but also we extended phase space. 

It is worthwhile to make some discussion regarding mass of the black hole in this case. Usually, the interpretation 
of mass is internal energy. In this point of view mass is a conserved quantity which is representing only the total 
internal energy. But in case of the new interpretation, mass of the black hole is a combination of internal energy 
and pressure. In other words, not only the mass of the black holes determines the internal energy and the shape of 
the black hole, but also we are expecting it to have information regarding interaction of constituents of black holes 
which is known as pressure of the system. Although, one may state that due to natural properties of black holes, it is 
impossible to study the interaction of constituents of black holes, one must take the approach of string theory to the 
matter into consideration. In other words, in context of string theory some attempts were made to study microstates 
of the black holes and their interpretation which can be extended by considering this point of view. 

The volume of the black hole is determined by topological structure of metric. Therefore, one expects that the 
calculated value of this quantity and the topological structure of the metric be in agreement. This result was obtained 
in calculation of volume. 

In linear electromagnetic field, the critical behavior of the system for both gravities were usual ones. One critical 
horizon radius was found, hence one phase transition was expected. The formation of the characteristic swallow tail 
was observed in G — T diagrams. In T — r-|_ and P — r+ diagrams the properties of phase transition were seen. But 
amazingly, in consideration of additional term the behavior of the system differed completely. 

By considering this nonlinear electromagnetic field, the structure of the phase diagrams and the thermodynamic 
behavior of the system were highly modified. Calculations regarding critical horizon radius lead to existence of 
two critical horizon radii. It is not unusual to find two critical horizon radii for black holes m , but in our 
case the presence of the second critical horizon radius was observed in plotted graphs. In other words, in case of 
this nonlinear electromagnetic field, contrary to other cases of nonlinear theories, a phase transition and a turning 
point were observed. The existence of these two points is related to existence of number of critical horizon radius. 
This result is evident by comparing Maxwell theory with this nonlinear theory. In case of Maxwell theory only one 
positive critical horizon radius was found which resulted into existence of one phase transition whereas for this case 
of nonlinear electromagnetic field, two critical horizon radii were found which resulted into existence of the phase 
transition and turning point. 

In case of G — T , three points were observed in which the behavior of the system changed in them. These points 
are representing different phases. In P — r-|_, for T = Tc, three horizon radii were found with same pressure (critical 
pressure) which indicates the existence of three phases and one turning point. For T — r+, in case of P = Pc, 
two horizon radii were found with same temperature (critical temperature). One of these horizons was located on 
subcritical isobar which is the usual phase transition point and the other one was out of the subcritical isobar. 

Another interesting issue was the effects of nonlinearity and GB parameters were opposite of each other. In case 
of nonlinearity parameter, the smaller critical horizon radius was an increasing function of nonlinearity parameter 
whereas the larger critical horizon radius was a decreasing function of /3. In this case the smaller critical horizon 
radius was highly function of nonlinearity parameter comparing to larger critical horizon radius. Therefore, in case 
of increasing nonlinearity parameter, the distance between these two critical horizon radii decreased. The critical 
temperature and pressure were increasing functions of nonlinearity parameter whereas Pcfc/Tc was a decreasing 
function of it. Interestingly, in case of GB parameter, the smaller critical horizon radius was a decreasing function 
of GB parameter and the larger critical horizon radius was an increasing function of it. The value of larger critical 
horizon radius was highly function of variation of GB parameter. Therefore, the distance between two critical horizon 
radii was an increasing function of GB parameter. As for the critical temperature and pressure, they were decreasing 
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functions of GB parameter. In this case, PcTc/Tc was also a decreasing function of GB parameter. 

These two parameters (a and /3) are describing two aspects of the black holes; gravitational and matter fields. 
In case of increasing /3, the nonlinearity behavior of the system increases. In context of electromagnetic tensor, 
increasing this parameter leads to decreasing value of this tensor ([5]). Therefore, the values of the electromagnetic 
tensor are decreasing functions of nonlinearity parameter. This fact leads to a conclusion that increasing the power 
of nonlinearity leads to decreasing power of electromagnetic field. But in case of GB parameter, by increasing its 
value, the power of the gravitational field increases. In case of this generalization due to considering higher orders of 
curvature scalar, we are increasing the gravitational force. As one can see, increasing gravitational force causes the 
critical temperature and pressure decrease which indicating that in this case phase transition is taking place in lower 
temperature, hence lower energy. Therefore, increasing the power of the gravitational force causes the system to have 
higher value of internal energy. This can also be seen in studying entropy of the GB gravity. As one can see, in case 
of spherical symmetric, entropy is an increasing function of GB parameter. In comparing GB and Einstein gravities, 
it is evident that GB gravity has higher entropy. Having higher entropy means that system has higher internal energy 
which causes the system to have phase transition in lower temperature. This result was found in studying phase 
diagrams of the GB and Einstein diagrams. 

By what was mentioned in last paragraph one can conclude following results: first of all generalization of GB gravity 
causes the internal energy of the system increases and the black hole has phase transition faster comparing to Einstein 
gravity. Therefore, in this theory black holes need to absorb less mass to have phase transition. Second, the entropy 
may be a function of the complexity of the black holes structure. Gonsidering higher orders of curvature scalar causes 
the system to have more complicated structure. Therefore, it is arguable that the entropy and the complexity of 
the structure of the system are related to each other. In other words, complexity of the structure is a measurement 
for entropy of the system. Nonlinearity of the system and gravitational force are two opposing factors. They are 
decreasing each others effects and at some points they may cancel each others effects. This argument is stating that it 
may be possible to fix parameters in a way which system has critical values same as ones in Einstein-Maxwell theory. 

Due to existence of the abnormal behavior in phase diagrams a question regarding the existence of reentrant phase 
transitions could rise which was answered through two methods. The coexistence diagrams and a new method showed 
that only one of the critical horizons are representing second order phase transition, while the other one could be 
another type of the phase transition. 

As for the effects of the dimensions, the smaller critical horizon, pressure, temperature and PcTc/Tq were increasing 
functions of dimensions whereas, the subcritical isobars and larger critical horizon radius were decreasing functions 
of it. This behavior is indicating that the higher dimensional black holes need to absorb more mass to have phase 
transition. It is worthwhile to mention that case n = 3 for Einstein gravity and n = 4 for GB gravity are considered 
to be special cases in these gravities. As for the GB gravity, in this specific dimension, some terms in Eqs. (TTSl) and 
(051) vanish. In case of Einstein gravity, it is evident from Eq. (1^51) that for case of n = 3 the power of the horizon 
radius is different from other dimensions in a way which causes the Gibbs free energy to have abnormal behavior 
which was observed in studying phase diagrams. 

As it was pointed out, considering whether we are working in context of stringy correction or classical Born-Infeld 
theories, the results may be different. As it was seen, in case of later correction, there was a limiting point which 
for specific values makes the system never acquire non-Van der Waals behavior, whereas in case of stringy corrected 
electromagnetic field such limitation did not exist and the system would have non-Van der Waals like behavior. This 
points out that in stringy corrected version of the theory, different terms have different dominant regions which 
correspondingly modify phase structure and phase transitions of the solutions. 

Finally, it is worthwhile to mention that regarding string theory, the quadratic Maxwell invariant arises independent 
of Born-Infeld electromagnetic field. If one conducts the study that was done in this paper by motivations of string 
theory (not correction to Maxwell theory), there will be no limitation for considered values. Therefore, one can omit 
vertical limiting lines in plotted graphs. In addition, we should note that modification in phase diagrams, previously 
was obtained by changing the gravitational part of action. Here in this paper, the modification of matter field caused 
the abnormal behavior of the phase diagrams. 
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